Exact Solution of selfconsistent Vlasov equation 
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An analytical solution of the selfconsistent Vlasov equation is presented. The time evolution is 
entirely determined by the initial distribution function. The largest Lyapunov exponent is calculated 
analytically. For special parameters of the model potential positive Lyapunov exponent is possible. 
This model may serve as a check for numerical codes solving selfconsistent Vlasov equations. The 
here presented method is also applicable for any system with analytical solution of the Hamilton 
equation for the formfactor of the potential. 



The selfconsistent Vlasov equation is one of the most 
frequently used equations for the time dependent de- 
scription of many-particle systems. Especially in nuclear 
physics this equation has been employed to describe mul- 
tifragmentation phenomena and collective oscillations. 
The numerical demands are appreciable to solve this 
equation in six phase space dimensions. It is apparently 
not widely known that there exists an analytical solv- 
able model from which the effects of selfconsistency can 
be studied. Here such a model is presented which shows 
that selfconsistency can lead to positive Lyapunov expo- 
nents. The explicit analytical solution provides a tool for 
checking numerical codes. 

The model single -particle Hamiltonian reads 



P 2 



2m 



(1) 



where V(r, t) is the mean field potential associated with 
the separable multipole-multipole force V1234 = ^312534 
resulting in 



V(P,t) = fig(f)Q(t). 



(2) 



This model has been employed e.g in Q| for numerical 
study of intrinsic chaoticity. We show in the following 
that this is indeed a consequence of selfconsistency. The 
selfconsistent solution requires 



Q(t) 



drdp 

(27T) 5 



9(r)f(r,P,t) 



(3) 



where the one-particle distribution function obeys the 
quasiclassical Vlasov equation 



dtf 



P 



d r f - d r Vd p f = 0. 



(4) 



I. METHOD OF SOLUTION 
A. Nonselfconsistent solution 

First we solve the differential equation (^) in nonself- 
consistent manner. It means we consider the time de- 
pendence of V due to selfconsistency as an external time 



dependence. Then the differential equation is a linear 
partial one and can be solved easily. 

We solve this equation by examining the differential 
equations for the equipotential lines. This can be found 
by rewriting (0) in the form of a sevendimensional gradi- 
ent 



(l,^,-d r V).df = 
m 



(5) 



where df = (f,d r ,d p ). Because of the fact that any 
gradient is perpendicular to the hypersphere we can 
see that any curve in this hyperplane, which can be 
characterized by a parametric representation in the way 
(p(s), r(s),t(s)), obeys the relation 



1 



Pi 
m 



-d n V = i(s) : Ti(s) : p\(s)\ 



x,y,z. 



(6) 



From this we can read off the differential equations for 
the hyperplane. Usually one eliminates the parameter s 
choosing the time t as a parameter. The result is the well 
known Hamilton equations 



d t p = -d r V(r,t) 



n -> r 

o t r = —. 

m 



(7) 



In the case we can solve this equations we would ob- 
tain a sixdimensional parametric solution of the Vlasov 
equation Cj = c.i(p,r,t); i = 1 — 6. Here the Cj are the 
integration constants of (Q). The general solution of the 
differential equation (^) is given as any function of these 
Ci. This function itself is determined by the inital distri- 
bution fo(po,fo). We reformulate the latter distribution 
therefore as a function of the Ci[po, r*o, t = 0] 



/o(Po, r Q ) = fc(c, t [p Q , r ,t = 0]) 



(8) 



which represents a variable substitution from (po,ro) co- 
ordinates into the new set of variables (ci). Therefore the 
initial distribution fo{pa,ro) is changed into f c (ci). The 
general solution of the Vlasov equation at any time can 
then be represented by 



f{p,r,t) = f c (ci\p,r,i\). 



(9) 
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We like to point out that instead of choosing the time as 
a parameter we have also the possibility to eliminate s 
in (H) by any variable p t or r, . This is especially helpfull 
for other models because then the energy appears as an 
explicit integral of motion. 



B. Selfconsistency 

Provided we know the nonselfconsistent solution (||) of 
the Vlasov equation ([!]) we can easily built in selfcon- 
sistency by employing (||). Introducing (||) into (p|) we 
obtain 



Q(t) = / drodpog(r[ci[r ,p }})fo{p ,r ) 



d(rp) 



d( Ci ) 



d{roPo) 
(10) 



which produces a complicated equation for Q and mo- 
ments of the initial distribution. If this can be solved, 
cq. (^) is the selfconsistent solution when Q(t) is intro- 
duced. 



II. ANALYTICAL MODEL 

Here we like to demonstrate the application of the 
method by an exactly solvable model. We choose a form 
factor of the form 



g(r) = a x x + a y y + a z z. 



For such a model system we can solve the Hamilton equa- 
tions exactly. This is performed by differentiating the 
second equation of (fj]) and inserting the first one 



d 2 ri aifi 



dt 2 m 
This is easily solved as 



Q(t); i = x,y,z. 



(12) 



= - a ^Q 2 {t) + c\t + <? i 



m 

m m 



(13) 



with 



Q 2 (t) = dt' dt"Q(t") 



such that the general solution of the Vlasov equation is 
any function of these c;. Taking the initial distribution 
as fo(p, r) into account corresponding to (||) we see that 
the general solution can be represented as 

f(p, r, t) = f (pi + fia l Q 1 (t), 

n H Q2(t) 

m 



-t(^ + ^Qi(t))). 

m m ' 



(16) 



One easily convince oneselves that this solution solves the 
Vlasov equation (||). 

The nonselfconsistent solution is given by writing the 
time dependence of Q(t) and reads 



/non(p, r, t) = f (pi + (XdtQt 

V a i ^.2 .Pi 



2m m 



(17) 



Now we employ the selfconsistency condition ([^J) in order 
to calculate Q(t). Therefore we shift coordinates in the 
integral to obtain the form 



Q(t) 



dpdr fi ai ,.,,Pi.s. ^ -« 
vai(ri Q 2 it) H t)j (p,r) 



a>i < n > Q2{t) H 1. 



(18) 



Double occurring indices i are summed over. Here we 
have introduced < a >o— J dpdfaf^p^r) / (2irh) 3 and 
the density n =< 1 > of the initial distribution fa. This 
selfconsistency condition ( |l8|) is solved by rewriting it as 
a differential equation 

Q"{t) = \ 2 Q(t) 
0(0) = at < n > =< r > 

Q (0) = di =< — > 

m m 



A = 



-an „ 
af 



rn 



(19) 



where the averaging <>o is performed about the initial 
distribution. The solution reads then 

v 1 

Q{t) = < r > coshAi+ < — > -sinliAi (20) 

m A 



Qt(t) = J dt'Q(t'). 
Rearranging now for Ci (p, r, t) we obtain 

4®f,t) = £* + ^ Ql ( t ) 

m m 

c 2 {p,r,t)= n + ^Q2{t) 
m 

-(- + —Qi(W 

m m 



(14) 



(15) 



from which one finds Qi, Q2 via (Q) and the selfconsis- 
tency solution ( |l6| ) follows. The selfconsistency solution 
is entirely determined by the initial distribution function. 
The further evolution is then given according to this ex- 
plicit time dependence. 

We see that we obtain in the case of n > 0, which 
means effective repulsive force, an oscillatory solutions. 
There is no chaotic behaviour. 
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The interesting solution is given by fj, < 0. There we 
have an exponentially decreasing — \J — and increas- 



ing mode y — ■ The later one defines indeed the 
largest Lyapunov exponent which is found to be positive 
here. This can be seen as follows. 

The mean momentum and position in any direction at 
a time t takes the form 



< Pi >t 

< n > t 



<L>o _^n Qi(t) 
m m 

<n>a + <^> t-^Q 2 (t). (21) 
m m 



Then we can calculate easily the mean phase - space dis- 
tance to an initial point via d = \J Ar 2 + Ap 2 from which 
we deduce the largest Lyapunov exponent A as 



For an explicitely solvable model the method is demon- 
strated. The largest Lyapunov exponent is calculated an- 
alytically and the conditions are investigated for the oc- 
currence of positive Lyapunov exponent. It is found that 
positive Lyapunov exponents are generated under certain 
potential parameters by selfconsistency. The model may 
serve as a useful check of numerical codes. 
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lim -In 



< p >f + < r >\ 



2 _ 



m 



M < 0. 



(22) 



With this expression we have presented a model which 
can be exactly solved within selfconsistent Vlasov equa- 
tion and shows explicitely that positive Lyapunov expo- 
nents are created by selfconsistency. 
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III. SUMMARY 



A method is presented to solve the selfconsistent 
Vlasov equation. The following recipie is proposed which 
is applicable for a model meanfield V — fig(r)Q(t) if the 
Hamiltion equations for this form factor g(r) arc inte- 
grable. 

1. Solution of the differential equations for equipoten- 
tial lines as a function of the nonselfconsistent (time 
dependent) potential. The solution is a parametric 
representation of the general solution. 

2. The initial distribution has to be expressed into 
this paramters. Then the time evolution is entirely 
determined by this parametric form of the initial 
distribution replacing the parameters by their time 
dependent form as derived in 1. 

3. The selfconsistency condition leads now to a gener- 
ally highly involved equation for the selfconsistent 
potential. This equation is derived using the non- 
selfconsistent solution of 2., which is a function of 
the potential by itself. 

4. Reintroducing this selfconsistent potential into the 
solution 2 the time evolution of the selfconsistent 
Vlasov equation is determined completely by the 
initial distribution. 
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